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ABSTRACT 

Observing the excitation mechanisms of incompressible transverse waves is vital for determining how energy propa- 
gates through the lower solar atmosphere. We aim to show the connection between convectively driven photospheric 
flows and incompressible chromospheric waves. The observations presented here show the propagation of incompress- 
ible motion through the quiet lower solar atmosphere, from the photosphere to the chromosphere. We determine 
photospheric flow vectors to search for signatures of vortex motion and compare results to photospheric flows present 
in convective simulations. Further, we search for the chromospheric response to vortex motions. Evidence is presented 
that suggests incompressible waves can be excited by the vortex motions of a strong magnetic flux concentration in 
the photosphere. A chromospheric counterpart to the photospheric vortex motion is also observed, presenting itself 
as a quasi-periodic torsional motion. Fine-scale, fibril structures that emanate from the chromospheric counterpart 
support transverse waves that are driven by the observed torsional motion. A new technique for obtaining details of 
transverse waves from time-distance diagrams is presented and the properties of transverse waves (e.g., amplitudes 
and periods) excited by the chromospheric torsional motion are measured. 

Subject headings: Sun: Photosphere, Sun: Chromosphere, Waves, MHD 



1. INTRODUCTION 

Magnetohydro dynamic (MHD) wave phenomena have 
now been observed to be ubiquitous throughout the solar 
atmosphere and are considered to be a potential mecha- 
nism for the transport of energy for the heating of the so- 
lar atmosphere an d the accelerat i on of t he solar wi n d (for 
reviews , see, e.g., lAschwandei^ . I2QQ5I : iKlimchukl . I2QQ6I : 
lErdelvi l2QQ6l ). However, no individual mechanism has 
yet been identified for converting a portion of the me- 
chanical energy generated in the Sun's convection zone 
to heat. 

Both incompressible and compressible MHD 
waves are widely reported through out the sola r 
atmo sp here (for reviews, se e, e.g.. iBaneriee et al 
2QQ7I. iZa aarashvih fc E rdelvii [2QQ1 iWand [201] 
De Moort el & Nakariakovi l2Ql"2[ ) and incomp r essibl e 
waves in the solar wind (e.g., iTu fc Marsch. 1995). 
Recent advances in ground and space-based solar 
telescopes have allowed for the detection of ubiquitous 
incompress ible transverse waves, both in t he chr o- 
mosphere (iDe Pontieu et all, I2QQ7I: iHe et al.L l2QQ9allbl : 



[19971 ). 

The solar magnetic field acts as a channel with the 
potential to distribute wave energy around the solar 
atmosphere. In the quiet Sun at the photospheric 
level, a large percentage of the magnetic flux appears 
to be concentrated into intense magnetic elements (spa- 
tial scales of 100 — 200 km) that outline the super- 
granule network. Some of these concentrations of mag- 
netic flux are observed to undergo a significant expan- 
sion at the chromosphere/Transition Region, forming 
magnetic funnels that could be preferential sites for so- 
lar wind acceleration or the legs of large- scale coronal 
loops (Peter, 2001). The remainder of the flux may 
form cell-spanning chromospheric struct ures, providing 
a ma gnetic canopy (Dowdv^t al., 1^8 6i i Rutteij , [2006L 
[20071 : IWedemever-Bohm et al.[ r [2009[ ). Smaller concen- 
tratio ns of magnetic flux al s o exist in the internetwork 
(e.g., Faurobert et al.[. r200l[: jDommguez Cerdeha et al.[ . 



Okamoto fc De Pontieu, '2011"; "AntoHn fc Verwichte", 
201 1[ : [kuridze et al, 2012; Morton e^al., 2012) 
and in the solar corona (|Tomczvk et all [2007L 
[Erdelvi fc Tarovanl [2008[ ). However, the incompressible 
waves are difficult to dissipate without strong gradients 
in A lfven speed (e.g., resonant absorption - ilonsorJ . 
[1978[ ) or some process to cause a cascade of the wave 
e nergy to higher frequencies (e.g., MHD turbulence 
- [Matthaeus et al.[ . [l999) where the y are efficientl 
dissipated by, e.g., cyclotron damping (|McKenzie et al 

Irichard . morton@northumbria. ac . ukl 



[msi; [Schrij ver fc Title, "2003') , potentially generating ad- 
ditional magnetic structure below the magnetic canopy. 

The largest source of wave-energy is the continu- 
ally evolving sub-photospheric convection, which gen- 
erates a wide spectrum of MHD fluctuations. The 
influence of convection on the magnetic elements is 
well studied at the photospheric level. Typically, 
the observations have focused on tracing the mo- 
tion of strong, localised magnetic flux elements, re- 
ferred to as magnetic bright points (MBPs) when seen 
in the G-band 4305 A line (e . g., [T itle et al. [. [1989 . 
Berge r et al., 1995, Ber ^erjt aLUl998r Berger fc Title . 
[20011 [Sanchez Almeida et al., 20M) or mo re recently in 
the H a line wings (Le ena arts et_aL, [2006L [Chitta etld] . 
[2012[ ). The continuous jostling of the MBPs leads to 
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the excitation of MHD waves and oscillations, well de- 
scribed by the t heory o f MHD oscillations in ma gnetic 
flux tubes (e.g., LS pruit' '19821 lEdwin fc Robert si. 119831: 
lErdelvi fc M"QrtQn . 2009; ErdelyL&Jedun, 2010^ 

One particular aspect of convective motion that has 
received increased attention recently is the generati on of 
vorti c es at the solar surfac e . Observations (i Brandt et al.l . 
119881 : iBonet et al.l . |2QQ8|; iBalmaceda et al... .2010.) and 
advanced simulations ( iMoll et al 1 [20TTI : IShelvag et al. l. 
l2011b ^ of solar granulation have revealed these vor- 
tex fea tures are almost everywhere on the solar disk. 



Wedemever-Bohm fc Rouppe van der VoortI (j2009l) and 
Wedemever-Bohm et al.l ()2012f ) have demonstrated that 



the vortex motion has an observable influence on the 
upper layers of the solar atmosphere, suggesting they 
may play a role in heating. Further, simulations of 
iFedun et al.l (|2011b) have demonstrated that the twist- 
ing of an open flux-tube by vortices at the photospheric 
level leads to the generation of MHD waves, in par- 
ticular incompressible (i.e., torsional Alfven and trans- 
verse/kink) waves. 

The transverse waves have recently been the subject 
of some controversy, with discussions on their nature, 
properties and nomenclature. The properties of the 
transverse waves in solar magnetic flux tubes are well 
understood (e.g. Edwin fc Rober tj, [jlgl lS2ruitl . 119821 : 
lEdwin fc Robertsl . I1983l lErdelvi fc Fedunl . |2007^, with a 
significant volume of research dedicated to developing 
realistic models, including, e.g., magnetic structuring 
(fVerth fc Erdelvil .2008: Ru der man et al.. 2008). densitv 



struc turing (jAndries et al.l . 120051: iDymova fc Ruderman 
2005 ), flux tube geometry fVan Doorsselaere et al 



2004 : pvmova fc Ruderman, 2006; Erdelyi fc Morton, 



J009, Morton fc Erdelvi, 200 9| ), multiple flux tubes 
([Robertson et al., 2010; Lun a et aJI. I2O IO Pascoe et aL 
201 iD . resonant da mping (iRuderman fc Roberts . 



"Goossens et all, I2002L iTerradas et al.l. 12006 



partial i onization (jSoler et al 



Terradas et al., 2010 . 

2009aiib : .Zaqarashvili et al.l . 1201 ID . 

One of the most important properties of the low- 
frequency tra nsverse waves is tha t they are highly in- 
compressible (iGoossens et al.l . 12009. ) , requiring large gra- 
dients in the Alfven speed to damp the waves. It is 
this property (shared with low-frequency bulk Alfven 
waves of a homogeneous plasma) that means incom- 
pressible waves can pass through the lower solar at- 
mosphere relatively undamped, making them a good 
candidate for acceleration of the solar wind. Obser- 
vations of incompressible waves in the quiet c hromo - 
sphere (jDe Pontieu et al.l. I2007L iMorton et al.l. I2012D 
quiet corona (jTomczyk et al L'2007', "Erdelvi fc Taroxad, 
I2OO8) and the solar wind (|Tu fc Marsch, 1995) support 
this idea. Another important property, that makes these 
waves relatively easy to observe, is that the transverse 
wave causes a physical displacement of the flux tubes 
axis. 

To assess exactly what role the incompressible 
MHD waves play in determining the dynamics 
of the solar atmosphere, one needs to combine 
advance models of wa ve propagation throughout 
the at mosphere (e.g ., iCranmer fc van Ballegooiien'. 
20051 ICranmer etaLl. Il o p, 



iM atsumoto fc Shibata, 
2010; "Antohn fc Shibata", '2010", Fedun et al., 2011a, 
van Ballegooijen et al.,. .2011. .Vigeesh et al.. . .2012.) with 



observations. The observations should aim to determine 
how the waves are being generated and track the wave- 
energy as it propagates through the solar atmosphere. 
This should reveal how much wave energy is present in 
each distinct region of the solar atmosphere and where 
the dissipation of this energy occurs. The tracking of 
waves through the solar atmosphere has been made a 
realistic possibility by the development of high spatial 
and temporal resolution multi-fllter systems, that can ob- 
serve different heights in the lower atmosphere simulta- 
neously (e.g.. Rapid Oscillations in the Solar Atmosphere 
(ROSA), CRisp Imaging Spectro-Polarimeter). The po- 
tential of these obser vational set-up s for w ave track- 
i ng is d emonstrated in IVecchio et al.l ()2007[ ) , iJess et ahl 
(20 12 a""!?). Such ground-based telescopes complement 
the satellites, e.g. Hinode, Solar Dynamics Observatory 
(SDO), which provide detailed observations in UV/EUV. 

This paper reports the observation of both chromo- 
spheric torsional Alfven and kink waves and provides 
evidence that the waves are generated by photospheric 
vortices. The analysed data is a high resolution, high ca- 
dence series of G-band, Calcium K and Hydrogen a lines 
observed with ROSA located at the Dunn Solar Tele- 
scope, USA. The ROSA observations are supplemented 
with SDO data. Photospheric flows are determined using 
Local Correlation Tracking (LCT) and compared to nu- 
merical simulations of convection (Section 4). A statisti- 
cal method for obtaining information on chromospheric 
kink waves is described based on the analysis of time- 
distance diagrams. 

2. OBSERVATIONS AND DATA REDUCTION 

The data presente d he re has been s t udied in part in 
iKuridze et all (j2012f ) and lMorton et al.l (j2012f ). The data 
were obtained at 15:41-16:51 UT on 29 September 2010, 
with the Dunn Solar Telescope at Sacramento Peak, 
USA. A six-camera system called ROSA was employed, 
details of which are given in Jess et al. (2010). In brief, a 
69''. 3 X 69'M region of the quiet solar atmosphere, posi- 
tioned close to disk centre (NO. 9, W6.8), was imaged with 
a spatial sampling of 0''.069 pixel"^. We note here this 
is the corrected pointing to that given in IKuridze et al.l . 
I20T2I . 

During the o bservations, high-order adaptive optics 
(jRimmeld . 12004) were used to correct for wave- front de- 
formations in real time. The seeing conditions were good 
but variable. Some frames in each of the data series are 
subject to signiflcant distortions. We exclude periods of 
seeing containing frames of inadequate quality. 

ROSA obtained images in multiple wavelengths includ- 
ing G-band (4305.5 A- 9 A width), Ca K (3933.7 A) 
narrowband (1 A) and Ha core (6562.8 A) narrowband 
(0.25 A) fllters. The G-band data was sampled at 16.6 
frame s~^. Ha at 2.075 frame s~^ and Ca K 1.66 frame 
s~^ and the image s were i mproved by using speckle re- 
construction ( Wog er" 3t al.l . I2OO8) utilising a 16 — 1 ratio. 
The cadence of the reconstructed G-band, Ha and Ca 
K time-series are 0.96 s, 7.7 s and 9.6 s, respectively. 
To ensure accurate co- alignment in all bandpasses, the 
broadband ti mes series were Fourier co-registered and 
de-stretched (|Jess et all . I2007D . 

Images from the SDO Atmospheric Imaging Assem- 
bly (AIA) are also presented to provide context for the 
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Fig. 1. — A region of the quiet solar atmosphere as observed by ROSA, (a) G-band image showing a 15 by 23 Mm^ sub-region of the 
ROSA field of view. Collections of magnetic bright points are clearly visible in the inter-granule lanes, (h) Ca K image of the lower 
chromosphere, (c) Ha core image of the mid to upper chromosphere. The existence of fine-scale structuring in the bandpass is evident, 
with both spicules/mottles and cell-spanning fibrils identifiable. The blue dots show the path along which perpendicular cross-cuts where 
taken, see Section 5 for further details. The bottom row displays images of the upper solar atmosphere observed by SDO/AIA. The panels 
show (d) 304 A, (e) 171 A and (f) 193 A filters. The axis show solar coordinates. The white boxes highlight the region which ROSA 
observed. 

state of the upper atmosphere above the above the ROSA 
field of view. The images from SDO are processed with 
aia-prep and have a spatial resolution of C)".6 pixel"^. 



3. THE MAGNETIC STRUCTURE 

The region of the Sun observed by ROSA is a typical 
quiet Sun area. The focus of the investigation will be on 
a sub-region of the ROSA field of view that contains a 
feature of particular interest. Images of this sub-region 
are shown in Fig. [TJ In the G-band images, a number 
of MBPs are clearly identifiable. The bright points are 
a good proxy for strong concentrations of magnetic fiux, 
however, strong magneti c fiux can still be pre sent with- 
out an associated MBP (jBerger fc Titld 1200 Specific 
conditions appear to be needed f or an MBP t o appear 
in G-band (e.g., .Carlsson et al.L |2004 llshikawa et al.. . 



|2007[ ). On viewing the movie of the G-band time se- 
ries, the MBPs are seen to be pushed and jostled by the 
convective motions of the granules, merging and splitting 
and appearing and disappearing from the bandpass. 

The Ca K images (Fig.[T|3) also show bright points, dis- 
playing a close mapping to the MBP's seen in G-band. 
This suggests the Ca K bright points are the chromo- 
spheric counterpart to the G-band MBPs. However, the 
Ca K bright points are more diffuse, possibly due to ex- 
pansion of the magnetic fiux with height and partially 
due to the properties of the filter. The bright points are 
also more persistent than the G-band MBPs, hence show 
a different morphology. The images also contain acoustic 
grains and reversed granulation, which have been iden- 
tified previously. These features (possibly including the 
diffuse MBP nature) are present due to the large width of 
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Fig. 2. — The left panel displays the G-band region enclosed in the red box in Fig. ^ The motion of an MBP flows derived with a 
point and click method is highlighted by the over-plotted stars (middle panel). The motion is from left to right. The right image is the 
synthesised G-band from the MURaM code. 

from SDO/AIA in the 304 A, 171 A and 193 A band- 
passes. A near-by active region is visible in the upper 
right hand corner of each image, however, it can be seen 
that the ROSA field of view is far from the active region 
in an area of quiet Sun. 



the filter bandpass used for these observations, meaning 
that the Ca K images here have contributions from both 
photospheric and chromospheric sourc es (for a detailed 
explanation, see Leenaarts et al.L [2006[ ). Selecting a suf- 
ficiently narrow filter to image Ca K line centre can avoid 
the photospheric contamination and allow mainly chro- 
mospheric con t ribut i ons (see, e.g., iCauzzT 
iReardon et all . 120091 : iPietarila et all . |2009). 

In the Ha core bandpass (Fig. [It), few signatures of the 
MBPs are visible. The region is covered with fine-scale 
structure, such as elongated, cell-spanning fibrils and a 
few shorter mottles (jet-like structures). Due to the nar- 
row width of the filter, the Ha core images show the 
'true ', magnetically dominated (low-/3) chromosphere 
^Ruttenl . 2006, 2007,). The observed features, i.e., the 
fibrils and mottles, are thought to highlight the chro- 
mospheric magnetic field, with the structures likely to 
be density enhancements along individual magnetic field 
lines (this is supported by recent m odelling of Ha line 
formation by Leenaarts et al.L |2012[ ) . 

In the following section we pay particular attention 
to the area highlighted by the dashed box in Fig. [TJ 
A number of fibrils and mottles can be seen to em- 
anate from this region. At the photospheric level, the 
G-band images reveal that a number of MBPs appear 
and disappear. The counterpart Ca K bright points re- 
main present, although some of the brighter elements 
fade and brighten with time. These underlying features 
lead us to suggest the highlighted Ha region is the chro- 
mospheric section of a of a magnetic fiux concentration, 
associated with an MBP, that has undergone significant 
expansion from the photosphere. The observed MBPs 
diameters are ~ 100 — 200 km while the Ha feature 
has a diameter of ^ 3000 km. The described magnetic 
structure is tho ught to be com mon in the quiet s olar at- 
mosphere (e.g., iriowdv et aO 119861 : iPeterl . I200l( ) and is 
also revealed in observations of chromospheric swirls by 
[Wedemever-Bohm & Rouppe van der Voort (2009) and 
[Wedemev er-Bohm et al. (2012). 

The magnetic structure above the chromosphere is 
much harder to determine. Fig. [T] (d-f) displays images 



4. FLUX TUBE DYNAMICS 

Now, our attention is focused on the dynamical be- 
haviour of the unique region identified in the previous 
section (dashed box Fig. [T]). To provide insight into the 
horizontal motions of the magnetic elements we exploit 
LCT. This technique was originally introduced into solar 
physics to track lar ge-scale (spatial and temporal) pho- 
tospheric fiows, e.g.. [November fc SimonI (p!988^) . and the 
technique has been developed further for use on magne- 
tograms by incorpora ting the induction equation (e.g., 
DAVE - [Scl^ . [2006[ ). 

We developed an LCT algorithm based on the rou- 
tine tr_get-disp.pro (developed by T. Tarbell) available 
in the SSW TRACE software tree. A recent version of 
this routine first apodizes a set of images using a Han- 
ning function, then computes the cross-correlation be- 
tween the images. The position of the maximum of the 
cross-correlation corresponds to the integer pixel shifts 
between the images. A local region centred on the max- 
imum is then fitted with a quadratic polynomial regres- 
sion equation to locate the maximum with sub-pixel ac- 
curacy, hence providing sub-pixel accuracy on the shifts. 

The routine was tested on sample data, which consists 
of hundreds of different images taken from the ROSA 
data series. The images were shifted using the routine 
shift -image. pro (available in the general SSW software 
tree). Taking averages over large numbers of images, we 
find that shifts as small as 0.03 pixels can be resolved to 
within ±0.005 pixels. The accuracy of the LCT improves 
with increasing shift size. There appears to be potential 
for great accuracy when calculating small shift values sta- 
tistically, however, a more rigorous test of tv-get-disp' s 
ability to resolve shifts is needed. 

From the derived shifts we are able to determine the 
horizontal photospheric fiow, Uh- Further, we can cal- 
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Fig. 3. — The top row shows the magnitudes of velocities of photospheric flows as determined from the observed G- band (left), the 
synthesised G-band (middle) and the MURaM code at a geometric height of 200 km above the continuum formation level (right). The 
velocities are averaged over 30 frames (~ 230 s). The middle row displays velocity vectors and divergence of photospheric flows. The 
arrows display the averaged velocity vectors determined with LCT from the observed (left) and simulated (middle) G-band. A histogram 
(right) showing the distribution of the horizontal velocity for the observed (red) and the simulated (black) data. The bottom row displays 
vorticity (s~^) calculated from the horizontal photospheric flows. The panels displays the averaged vorticity from the observed (left) and 
simulated (middle) G-band. A histogram (right) showing the distribution of the vorticity for the observed (red) and the simulated (black) 
data. 



6 



Morton et al. 






Fig. 4. — The chromospheric response to photospheric vortices. Torsional motion is observed in Ho; and demonstrated by applying LCT 
to frames between t = 207 — 277 s (left) and t = 610 — 690 s (right). The torsional motion is centred approximately on the point located at 
(3.8,4.5) Mm. The vectors are over plotted on the unsharp masked and atrous filtered Ha data at t = 200 s (left) and t = 600 s (right). The 
region shown corresponds to the dashed box in Fig. [T]:. The bottom row displays enlarged versions of the regions centred on (3.8,4.5) Mm 
(as marked by the boxes). The density of vectors has also been increased by a factor of two and a half. 
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Fig. 5. — Overview of the technique for identifying transverse waves. The left panel shows time-distance plot of a cross-cut taken 
perpendicular to the fibrils (i.e., perpendicular to the dotted blue line in Fig. ^). The middle panel is the same cross-cut that has been 
subject to unsharp mask and with the high frequency component removed. The right panel displays the fibril paths obtained from the 
time-distance plot in the middle panel. The lines highlight the dark features. 



culate the divergence of the velocity, V ■ Uh^ and the 
photospheric vorticity, V xUh- 

4.1. LCT on G-band data 

The G-band data series has a very high cadence at 
0.96 s. However, the previously reported values for pho- 
tospheric flows are 1 — 2 kms~^, hence the expected 
shift between each frame is < 0.03 pixels. To improve 
accuracy, every eighth frame is then selected from the 
G-band series, increasing the cadence to 7.68 s and ex- 
pected shifts to > 0.08 pixels (0.5 kms~^). Further, a 
sonic filter is applied to the data to suppress the influ- 
ence of p-modes and stochastic variations in intensity. 
The data is then re-sampled using linear interpolation to 
achieve a pixel size of 25 x 25 km^. 

The LCT is performed on frames chosen from the first 
sixty frames because the seeing during this period is rela- 
tively stable. For each data set that we apply LCT to, the 
LCT is preformed on many sub-regions of the overall im- 
age, with window sizes determined by the dominant fea- 
tures at that wavelength. A window size of 40 x 40 pixel^ 
(1x1 Mm^) is used such that features on the order of 
the granulation will contribute to the cross-correlation. 
This, along with the apodization of the window, will re- 
duce the effect of large intensity gradients on resolving 
shifts. 

The 7.5 X 7.5 Mm^ boxed region in Fig. [TJi is selected 
(close-up of this region is shown in Fig. [2] left panel) and 
subject to the LCT algorithm. The photospheric flow is 
calculated over a 30 frame average, revealing long-lived 
flow patterns. 

The magnitude of the photospheric velocities are 
shown in Fig. [H A histogram plot of the velocity distri- 



buti on (Fig;. [3]) is in line with results fr o m MB P tracking, 
e.g., " Kevs et all (j2011f ): iChitta et al7 (2012) , and from 
l ong-t er m flow meas urements, e.g., Novem ber fc Sim"Qnl 
(jl988f ). iBerger et aP (p98). The velocity vectors and 
divergence are also shown in Fig. [3l The directions of 
the velocity vectors results are consistent with the flow 
direction as seen by a visual inspection of the G-band 
movie. Evidence for vortices are present, with the veloc- 
ity vectors marking out curved paths or arcs thought to 
be demonstrating the vortex behaviour. This is verified 
by calculating the vorticity from the horizontal velocity, 
which demonstrates that vorticity is present across the 
the photospheric region (Fig. [3]- bottom row, left panel). 
A histogram of vorticity values is given in Fig. [3] and 
the values obtained here are in agre ement with previous 
observational measurements ( e.g;., iBrandt et al.L 119881 : 
I Vargas Dommguez et al.l . l2011f ). 

In addition, the motion of an individual MBP is 
tracked over twenty frames using a crude (but efficient) 
point and click method. The MBP selected lies directly 
under the chromospheric feature, which suggests it is the 
photospheric counterpart. The tracking of the individual 
bright point reveals the MBP follows an arc-shaped path 
across the photosphere, demonstrating the MBP is in the 
presence of a photospheric vortex. This is confirmed as 
the MBP is also seen to lie in a region of strong, negative 
vorticity in Fig. [3l 

4.2. LCT on simulated G-band data 

To validate the results of the LCT technique, we per- 
formed tests on a numerically simulated G-band image 
sequence. The synthetic data were produced using ra- 
diative diagnostics of the three-dimensional numerical 
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model of magnetised lower solar atmosphere (i.e. con- 
vection zone and p hotosphere) generat ed by MURaM ra- 
diative MHD code (IVogle r et al.Ll2QQ5[ ). The domain size 
chosen isl2xl2xl.4 Mm^ in the horizontal and verti- 
cal directions respectively, and the horizontal resolution 
is 25 km. The detailed G-band radiative diagnostics was 
carried out using the radiative transport code, designed 
for computationally efficient wide band filter calculations 
(see, Jess et al., 2012b for more details). We have gener- 
ated 800 G-band images of solar granulation with the 
average vertical magnetic flux of 200 G (the value of 
magnetic fiel d strength is consistent with recent obser- 
vations, e.g., iQrozco Suarez fc Bellot Rubid [20121 ) . The 
time cadence for the image sequence is about 2 seconds 
and not constant due to variation of the time step in the 
numerical simulation. 

To provide a direct comparison with the observational 
analysis described in the previous section, the resolution 
of the simulated data is reduced to 50 km, using every 
fourth frame to achieve a similar cadence. We apply 
the sonic filter and re-sample using linear interpolation 
to 25 km. The LCT results are averaged, again, over 
30 frames. The results of LCT on the simulated G-band 
image sequence demonstrate qualitatively similar plasma 
flow patterns in both the granular and intergranular re- 
gions to the observed data (Fig.[3j). The results from the 
simulated data show slightly greater values of divergence 
and vorticity. The absolute values of velocity and vor- 
ticity in the simulated data are plotted in a histogram 
(Fig. [3]) for direct comparison between the observed and 
synthesised G-band images. The distributions can be 
seen to be comparable although slight differences exist 
between the observed and simulated values. This can be 
attributed to the fact that the observed and simulated re- 
gions are inherently different regions, with an unknown 
quantity of magnetic flux in the observed region. Some of 
the differences may also lie in the quality of the observed 
data. A measure of the difference in quality is the root 
mean squared contrast. The synthesised G-band images 
have a mean RMS of 0.17 ±0.003, while the observed G- 
band images have a much lower mean RMS of 0.11 ±0.02 
and much greater variation in RMS between frames. 

The flow speeds in the simulations have also been mea- 
sured at the geometrical height of approximately 200 km 
above the continuum formation level. They show the 
presence of the flow features similar to those detected by 
LCT applied to the simulated G-band images (Fig. [3j). 
However, due to measurement of the flow speeds at the 
constant geometrical height instead of at optical depth 
corresponding to the G-band radiation formation, only 
a qualitative comparison of these speeds with the speeds 
obtained by LCT is valid and possible. 

The simulations are known to be replete with vor- 
tex motions. The agreement seen between LCT results 
performed on the observed and simulated data gives 
us confidence that the photosphere also contains nu- 
merous vortex structures. Torsional Alfven waves have 
been shown to be excited at the photospheric level in 
a number of numerical simulations. They are excited 
in magnetic photospheric vortex structures by random 
motions, such as granular flows (see, for example, nu- 
mer ical simulations of photospheric magneto-convection , 
e.g. Ivan Ballegooiien et al.L[2011l : lKitiashvih et al.Ll2011l : 



Shelyag et all |2011al, i Shelvag et all . I2012al lbl: iMoll et al.L 
20121 ). and by syntheti c photospheric t or sional (vortex- 
type) drivers (e.g., Fed un et al.L l20l"lal l3. IVigeesh e t aTl . 
I2OI2I ) in the idealised simulations of solar magnetic flux 
tubes. 

We note the magnetic vortex structures significantly 
differ in their nature and dynamics from th e non- 
magne tic vortices (see , e.g., iStein fc Nordlundl . I1998L 
Kitias hvili et al.L l2012f ) . They, and the torsional Alfven 
waves they produce, generate large amounts of Poynt- 
ing flux (see, e.g., Fedun et al., 2011b, Shelv ag et al.L 
2012a.b). As has been shown by iWe demeyer- Bohm et al.l 
(2012), the magnetic vortex structures expand from the 
photospheric level into the corona. They are suggested 
channels for electromagnetic energy to tunnel into the 
higher layers of the solar atmosphere and, therefore, such 
structures are considered as potential channels of energy 
transport for coronal heating. 

4.3. LCT on Ha data 

On inspecting the Ha movie, a torsional (or rota- 
tional) motion is observable in the area highlighted 
with the dashed box in Fig. [ih. A chromospheric re- 
sponse to photospheric vortices has already been ob- 
served in high-resolution observatio ns of the lower chro- 
mosp here (Wedemever-Bohm & Ro uppe van der Voortl . 
I2009L IVMemeyer-Bohm et al., 2012^. So far, the re- 
ported chromospheric swirl events appear to demonstrate 
rotation in one direction only. Here, the Ha region ap- 
pears to ffist rotate in one direction and then rotate in 
opposite direction quasi-periodically. Such behaviour has 
been observed in the intergranul ar vortices in r ealistic 
simula tions of magneto-convection IShelvag et al.] ()2012aL 
[20TTbl ). as weh as in the idealised simulations of flux 
tubes in the quiet Sun (Fedun et al., 2011c). This would 
suggest the presence of torsional Alfven waves, which 
have also been identified pre viously in chromo spheric fea- 
tures situated above MBPs (|Jess erall . l2009l ). 

To demonstrate this behaviour, LCT is applied to the 
Ha data. Before performing LCT, the unsharp-mask 
(USM) technique is applied to the full field of view and 
the mask is subtracted. The mask is created using a box- 
car average of width 20 pixels. This technique has two 
benefits: (i) the removal of large spatial-scale intensity 
perturbations which could influence results; (ii) an in- 
creased contrast between chromospheric features and the 
backg round. The intensity perturbations (Morton~el^, 
are removed from the time series because they are 
coherent over a large spatial- scale. Further, an atrous 
spatial ffitering algorithm (e.g. JStarck fc Murtaghl . l2QQ2h 
is applied to each frame. The highest frequency com- 
ponent from the spatial filtering contains mainly noise, 
which is then subtracted to improve the signal to noise 
in each image. 

The torsional motion appears quasi-periodic with an 
apparent period between 120-180s, hence, the LCT is 
averaged over sets of 9 consecutive frames. Note, the pe- 
riod range given is based on the apparent rate at which 
the motion appears to change its direction of rotation. 
The time-series shows the rotation change directions a 
number of times over it's duration. This is calculated 
by eye from viewing the Ha time-series and we do not 
consider it a rigorous measurement of the period. The 
velocity vector plots in Fig. |4]conffim the visual impres- 
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sion of torsional motion. The first vector plot is averaged 
over the frames in the range 207 — 277 s and the second 
vector plot is averaged over 610 — 690 s. These time 
ranges are used as they correspond to the times when 
the torsional motion is most evident. The torsional mo- 
tion is centred on the region around (3.8,4.5) Mm in the 
plots and a zoom of this region (Fig. |4] - bottom row) 
reveals the signature of the torsional motion. The max- 
imum values of the velocity amplitude of the torsional 
motion are ^ 7 kms~^ with an average and standard 
deviation of 1.8 ± 1.4 kms~^. These values are essen- 
tially the time-averaged velocity amplitudes, while the 
peak amplitudes are likely greater as the peak value is a 
factor of a/2 greater than the time-averaged amplitude. 

5. TRANSVERSE WAVES 

The focus of this work is now shifted to the fine-scale 
structure that originates in the previously identified Ha 
region. The fine-scale structure we are interested in are 
the dark, absorption features, which are mainly elon- 
gated, inclined fibrils and a few mottles (nearly vertical 
structures thought to be similar to spicules). On inspect- 
ing the movie of Ha, the fine-scale structure can be iden- 
tified to exhibit transverse motion. The fibrils are rooted 
in the region where the torsional motion is observed and 
the torsional motion shakes the fibrils footpoints side to 
side, which drives the transverse waves. In this section we 
study the transverse waves supported by these structures 
and demonstrate new techniques for obtaining informa- 
tion about the waves. 

To model these transverse waves, the assumption is 
made that the waves are propagating along an over- 
dense, magnetic fiux tube embedded in an ambient 
plasma environment. The speed at which these waves 
propagate is known as the kink speed, which is defined 
as, 

c. = ,/S5±^. (1) 

V Pi + Pe 

where p is the density and va = 5^/(p/io) is the Alfven 
speed, B is the magnetic field and fio is the magnetic per- 
meability. Here, we follow convention labelling internal 
and ambient plasma parameters with subscripts i and e, 
respectively. 

The first step in the analysis of the observed trans- 
verse waves is to apply the USM procedure to the Ha 
images. The orientation of the fine-scale structure is de- 
termined (see dotted line over plotted on the Ha image. 
Fig. [It) and a series of 42 cross-cuts, separated by 2 
pixels, are placed perpendicular to their axis. From the 
cross-cuts, time-distance plots are created, an example 
of which is shown in Fig. [5l The time-distance plots 
clearly reveal the transverse motions of the fibril struc- 
tures. Next, to reduce the background noise, an atrous 
filtering algorithm is applied to the time-distance plots 
and the highest frequency component is subtracted. 

The individual, dark fibrils in each cross-cut are then 
located in the following manner. Cross-sectional flux pro- 
files of fibrils have an almost Gaussian shape, with the 
minimum value of intensity at the centre of the fibril 
cross-section. For each time slice, the pixels with the 
minimum intensity in a localised region are located. A 
gradient threshold is then applied allowing for isolated 



fibrils to be located. The gradient of the neighbouring 5 
points on either side of the minimum value pixel is cal- 
culated and has to be larger than the threshold value 
for the minimum point to be selected. Should two fibrils 
approach each other or begin to cross, the gradient be- 
comes shallower around the minima. Optimisation of the 
threshold level allows for a significant number of isolated 
fibrils to be identified (Fig. [Sb), with the fibril paths in 
time being traced out. This method only locates the cen- 
tre of the structure to within one pixel and hence has an 
error of ±50 km on each point. 

For each fibril path identified at this stage, we then 
select only the paths that obey certain conditions: (i) 
the path has to be longer than 5 pixels; (ii) The max- 
imum time separation between neighbouring points in 
the thread has to be less than 4 pixels (30 s), the 'miss- 
ing' pixels are then recorded as having the same spatial 
position as that of the first of the time-separated neigh- 
bours; (iii) Points neighbouring in time cannot be sepa- 
rated spatially by more than 4 pixels (200 km). This last 
condition would limit any measured velocity amplitudes 
to < 26 kms-^ 

Once we have selected the fibril paths, the results are 
then fitted with a Levenberg-Marquard t non -linear fit- 
ting algorithm (mpfit.pro - Markwardt. , |2009[ ). A func- 
tion of the form 

F{t) = G{t) ^ Asm{ujt - (^) (2) 

is used to fit the oscillations. Here, G{t) is a linear func- 
tion, A is the displacement amplitude of the oscillation, 
u the frequency and (j) the phase of the oscillation. The 
fitting algorithm is supplied with the errors on each data 
point, where it is assumed the given error is the one-a 
uncertainty. The fitting routine also calculates the one-cr 
error to each fit parameter. The fitted results are kept 
if the threads are longer than 3 • 27r/4a; and the errors 
on the fitted parameters are smaller than the magnitude 
of the parameter values. Transverse velocity amplitudes 
for the waves can be obtained from the fit using the re- 
lation v± = ujA and errors in the velocity amplitude are 
calculated by summing in quadrature. 

Over the 42 cross-cuts a total of 1100 oscillations are 
measured. A histogram of the periods, transverse am- 
plitudes and velocity amplitudes are given in Fig. [6l A 
Gaussian fit to the distributions gives means and stan- 
dard deviations of P = 94 ± 61 s, A = 71 ± 37 km and 
v± = 4.5 ±1.8 kms~^. It is also interesting to plot period 
versus displacement amplitude (Fig. [6]- bottom row left 
panel). It is apparent that there is a direct relationship 
between transverse displacement amplitude and period. 
We perform a weighted linear fit to the data points, cal- 
culated with the Levenberg-Marquardt algorithm with 
the constraint the line passes through (0,0). The result 
is over-plotted and the derived velocity amplitude from 
the trend is ~ 4.2 kms~^, in agreement with the mean 
value. 

In Fig. [6] (bottom row right panel) the velocity ampli- 
tude as a function of period is shown. Qualitatively the 
velocity amplitude appears relatively constant as a func- 
tion of period, however, a weighted linear fit to the data 
points suggest the velocity amplitude decreases with pe- 
riod. This implies that the energy, which is proportional 
to the square of the velocity, contained within this fre- 
quency range of waves (2-30 mHz) is approximately con- 
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Fig. 6. — Histograms of measured properties of transverse motions for chromospheric fine structure originating in the region displaying 
torsional motion (top row). The histrograms show, from left to right, the period, transverse displacement amplitude and transverse velocity 
amplitude. The mean and standard deviations are P = 94 ± 61 s, A = 71 ± 37 km and v± = 4.5 ib 1.8 kms"-*^. Measured properties of 
the transverse waves (bottom row), (left) Period versus displacement amplitude. The solid line is a weighted fit to data points that gives 
a velocity amplitude of ~ 4.5 kms"-*^. (right) Period versus velocity amplitude. The dashed line is weighted linear fit to the data points. 



stant as a function of period. It is worth noting that the 
estimates of the temporal power spectra of the horizon- 
tal velocity of the photosphere by Matsumoto & Kitai 
(12010) using the G-band filter of Hinode/SOT showed a 
change of trend in the log-log gradients of Fourier power 
vs frequency at 4.7 mHz (approximately 213 s period). 
For frequencies less than 4.7 mHz they found an f~^'^ re- 
lationship and for frequencies greater than 4.7 mHz it was 
y-2.4 rpj^^ derived photospheric power spectra, found 
by the method of LOT were taken as a proxy for the 
actual photospheric d river of transverse wave s for a nu- 
merical simulation by Matsu moto fc Shibatal ()2QlQf ) (in 
their case the m = torsional Alfven wave and not the 
m = 1 kink waves we observe here). A number of ques- 
tions now arise, in comparing our measurements of trans- 
verse chromospheric velocity amplitude vs period mea- 
surements in chromospher ic fibrils with the photospheric 
power spectra of Matsumo to fc Kit ai (2010). Firstly, our 
results suggest that Mats umoto fc Shibataj (|2Q1Q>) may 
have used the incorrect input power spectra, i.e., we see 
no evidence of wave energy decreasing with decrease in 
period. This is particularly relevant for wave heating 
since, mechanisms such as resonant absorption, phase 
mixing and ion-neutral damping are all more efficient 
at higher frequencies. On the other hand, if the photo- 
spheric power spectra obtained by Matsumoto fc Kitai 
([2010 ) are a reasonable proxy for the driver of transverse 
waves, the output is notably different at chromospheric 



heights and must be explained by further study. 

We should note that all fits calculated here are under 
the assumption that the error in the period can be ne- 
glected over the error in the transverse displacement and 
velocity amplitudes. This may be justified as the mean 
error on the period, transverse displacement and velocity 
amplitudes are 8%, 21% and 23%, respectively. 

6. DISCUSSION AND CONCLUSIONS 

Tracking wave propagation through the lower layers of 
the solar atmosphere is an important step in being able to 
assess the viability of the various suggested mechanisms 
for atmospheric heating and solar wind acceleration. In 
Section. 4 we present observational evidence that photo- 
spheric vortices occur in the regions of strong magnetic 
flux concentrations. This is found to be in good agree- 
ment with the predictions of advanced models of solar 
magnet-convection. Further, we provide evidence that 
demonstrates the photospheric vortices can excite tor- 
sional motions in the chromosphere. This is supported by 

previous results from v ar ious n u merical simula tions (e.g., 
■ ii rr ^ 



van Ballesooiien et al., '2011"; 'Kitiash vih et al.L 1201 ll : 
Shelvag et al., 2011a, 2012a.b, Fedu n et al.L I2011al l3. 
iMoh et al.L 1201 2[ lYTgeesh et al.L I20T2I) and previ- 
ous observations (e.g., IWedemever-Bohm et al.L |2012[ ). 
The observation of the vortex and torsional motions 
is achieved by exploiting LOT to track the mo- 
tions of both photospheric and chromospheric fea- 
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tures and identifying the signatures in the result- 
ing velocity vector plots. The chromospheric tor- 
sional motion appears periodic here (unlike the uni- 
directional motion identified in swirl e v ents, e.g., 
IWedemever-Boh m fc Rouppe van der Voorti . ^00§), al- 
though we cannot resolve the periodic behaviour explic- 
itly. The vector plots demonstrate that the region rotates 
one way and, at a later time, in the opposite direction. 
A better way to obtain details on torsional waves may 
be to exploit spectroscopic observations (e.g. JJess et al.l . 
^009). 

Further, the torsional motion is observed to excite 
transverse waves in the chromospheric fine-structure 
whose footpoints are rooted in the region. The fine- 
structure takes the form of elongated absorption fea- 
tures that are fibrils and a few shorter features, possi- 
bly mottles. The transverse velocity amplitudes of the 
transverse waves are comparable to the velocity ampli- 
tudes obtained from averaging of the torsional motion. 
The observations imply that the transverse waves are 
driven by the torsional motion, which in turn is driven 
by the photospheric vortices. The transfer of energy 
from photospheric vorticies to Alfven waves and chromo- 
spheric transverse waves has been demonstrated in nu- 
merical simulations of the lower solar atmosphere (e.g., 
iFedun et al.ll2Qllb[ ). 

The measured wave properties also reveal interesting 
information about the chromospheric waves. In partic- 
ular, the measured velocity amplitudes (Fig. [6]) suggest 
that waves with periods in the range 50 — 500 s have sim- 
ilar power. However, the results presented here are re- 
stricted to incompressible waves confined in one particu- 
lar structure in the chromosphere. A more in depth study 
taking into account numerous chromospheric structures 
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is required before more definite conclusions on this issue 
can be drawn. 

Finally, we highlight here that our results have some 
limitations. The limitations are related to the inability 
to resolve small amplitude (< 25 km) transverse waves 
due to spatial resolution. The demonstrated relationship 
between displacement amplitude and period (Fig. [5]) im- 
plies that spatial resolution also limits the ability to de- 
tect short period waves (< 30 s). On the other end of 
the scale, the known lifetimes of the chromospheric struc- 
tures (i.e., 3 — 5 minutes) also limits the observations of 
longer period waves (> 200 s). This is reflected in the 
low number of observations of wave with longer periods. 
Some of these limitations may be overcome if more ad- 
vanced techniques are employed for wave tracking, e.g., 
calculating sub-pixel displacements could reveal smaller 
amplitude wave phenomena. 
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